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FOREWORD 


Opinions,  interpretations,  conclusions  and  recommendations  are 

those  of  the  author  and  are  not  necessarily  endorsed  bv  the 
U.S.  Army. 


Where  copyrighted  material  is  quoted,  permission  has  been 


obtained  to  use  such  material. 

l^ere  material  from  documents  designated  for  limited 

distribution  is  quoted,  permission  has  been  obtained  to  use  the 
material. 

^^^^^^^^Tcitations  of  commercial  organizations  and  trade  names  in 
this  report  do  not  constitute  an  official  Department  of  Army 

endorsement  or  approval  of  the  products  or  services  of  these 
organizations. 

In  conducting  research  using  animals,  the  investigator (s) 
adhered  to  the  "Guide  for  the  Care  and  Use  of  Laboratory 
Animals , "  prepared  by  the  Committee  on  Care  and  Use  of  Laboratory 
Animals  of  the  Institute  of  Laboratory  Resources,  National 
Research  Council  (Nil,  Publication  No.  86-23,  Revised  1985). 

For  the  protection  of  human  subjects,  the  investigator (s) 


adhered  to  policies  of  applicable  Federal  Law  45  CFR  46 

— In  conducting  research  utilizing  recombinant  DNA  technology 
the  investigator (s)  adhered  to  current  guidelines  promulgated  by" 
the  National  Institutes  of  Health. 

In  the  conduct  of  research  utilizing  recombinant  DNA,  the 


investigator (s)  adhered  to  the  NIH  Guidelines  for  Research 
Involving  Recombinant  DNA  Molecules. 

In  the  conduct  of  research  involving  hazardous  organisms. 


the  investigator (s)  adhered  to  the  CDC-NIH  Guide  for  Biosafety  in 
Microbiological  and  Biomedical  Laboratories . 


13.  ABSTRACT: 


The  goal  of  our  proposed  work  was  to  improve  the  breast  tumor  diagnostic 
capabilities  of  optical  spectroscopy  and  imaging  based  on  diffusing  near  infrared 
light.  We  introduced  new  theoretical  methodologies  for  image  reconstuction 
within  the  human  breast,  and  planned  to  develop  these  theoretical  approaches 
further,  and  test  these  approaches  experimentally.  The  primary  benefit 
derivable  from  these  ideas  is  improved  breast  image  fidelity,  which  in  turn, 
would  enable  clinicians  to  fully  exploit  the  new  spectroscopic  and  scattering 
contrast  mechanisms  available  with  the  optical  method  for  increased  tumor 
sensitivity  and  specificity.  Since  submission  of  this  proposal  we  have 
experimentally  demonstrated  the  theoretical  ideas  underlying  the  near-field 
Fast  Fourier  Transform  approach.  These  experiments  produced  two-dimensional 
projection  images  of  thin  objects  embedded  in  infinite  media.  The  apparatus 
however,  was  not  well  suited  for  realistic  parallel  plate  soft  breast  compression 
geometries,  so  during  the  last  year  we  have  completed  the  construction  of  an 
instrument  that  is  far  better  suited  for  these  more  realistic  studies.  We  have 
also  explored  various  theoretical  issues  including  the  effects  of  imaging  filters  on 
the  reconstructions,  and  the  potential  for  three-dimensional  reconstruction. 


14.  Subject  Terms;  Photon-Migration,  Photon-Density- Waves,  Near-Infrared- 
Spectroscopy,  Fast-Fourier-Transform,  Image-Reconstruction,  Breast  Cancer 
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INTRODUCTION 


The  application  of  the  optical  method  to  breast  tumor  imaging  and 
specification  is  attractive  for  several  reasons.  The  techniques  utilize  non-ionizing 
radiation,  are  non-invasive,  and  are  often  technologically  simple  and  fast.  These 
features  have  lead,  in  some  cases,  to  compact,  portable,  and  inexpensive  clinical 
instruments  in  other  fields  of  medicine.  In  addition  the  optical  method  has  a 
great  range  of  measurable  parameters  that  could  enhance  tumor  sensitivity  and 
specificity.  For  example,  the  higher  metabolic  activity  of  the  rapidly  growing 
tumor  increases  organelle  population,  particularly  mitochondria,  which  leads  to 
an  increasing  scattering  factor  for  the  tumor.  Blood  dynamics  such  as  oxygen 
saturation  are  expected  to  be  substantially  different  in  the  rapidly  growing 
tumor,  and  will  idter  optical  absorption  factors;  the  optical  absorption, 
fluorescence,  and  scattering  of  contrast  agents  such  as  indocyanine  green  that 
occupy  vascular  and  extravasular  space,  may  also  provide  useful  forms  of 
sensitization. 

Optical  characterization  of  the  breast  has  been  attempted  since  1929 
[Cutler,  1929]  when  the  term  “diaphanography”  was  applied  to  shadow  graphs  of 
breast  tissue.  Although  the  method  has  been  tested  over  the  years  on  some 
2000  patients,  the  method  was  largely  inadequate  for  clinical  use  since  it  was 
difficult  to  separate  the  effects  of  absorption  and  scattering  within  the  tissue. 
Typical  photon  random  walk  steplengths  in  tissues  are  in  the  range  of  1  mm 
(corresponding  to  a  reduced  scattering  coefficient  oflO  cm'‘)  for  photons  in  the 
near  infrared  (i.e.  ~800  nm).  Even  in  the  region  of  low  tissue  absorption  (i.e. 
between  600  nm  and  1300  nm),  this  high  degree  of  scattering  tends  to  distort 
spectroscopic  information  and  blur  optical  images  as  a  result  of  the  large 
distribution  of  photon  pathways  through  the  tissue  that  contribute  to  measured 
signals. 

A  recent  and  critical  advance  in  the  photon  migration  field  however,  [see 
papers  and  references  in  Chance  &  Alfano,  1995;  Yodh  &  Chance,  1995],  has  been 
the  recognition  and  widespread  acceptance  that  light  transport  over  long 
distances  in  tissues  is  well  approximated  as  a  diffusive  process.  Using  this  basic 
physical  model,  it  is  possible  to  quantitatively  separate  tissue  scattering  effects 
from  tissue  absorption  effects,  and  one  can  accurately  incorporate  the  effects  of 
boundaries,  such  as  the  air-tissue  interface,  into  the  transport  theory  [Patterson 
et  al.,  1991a,  1991b].  Waves  of  diffuse  light  energy  density  [Gratton,  1990]  (or 
their  time-domain  analogs  [Patterson  et  al.,  1989;  Delpy  et  al,  1988;  Jacques, 

1989;  Benaron  &  Stevenson,  1993] )  are  generated,  and  propagate  deeply  in 
tissues  (i.e.  as  much  as  10  cm)  while  obeying  simple  optical  rules  such  as 
refraction  [O’Leary  et  al.,  1992],  diffraction  [Boas  et  al.,  1993;  Fishkin  &  Gratton, 
1993],  interference  [Schmitt  et  al.,  1992],  and  dispersion  [Tromberg  et  al.,  1993; 
Sevick  et  al.,  1992]  as  they  encounter  variations  in  tissue  optical  properties. 
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The  promise  of  the  optical  method  has  led  to  the  first  clinical  studies  of  the 
optical  technique  in  the  human  breast.  In  these  experiments  [Franceschini  et  al., 
1996;  Fantini  et  al.,  1996]  diffusive  waves  were  applied  to  the  human  breast  in 
the  transmission  mode  in  a  conventional  X-ray  “hard  compression”  parallel  plate 
source-detection  geometry.  The  team  obtained  wave  amplitude  and  phase-shift 
data,  and  then  processed  Ae  data  using  a  phase-corrected  linear  backprojection 
approach  [Fantini  et  al.,  1996].  The  inversion  procedure  did  not  separate 
scattering  from  absorption,  and  furthermore,  the  inversion  model  was  extremely 
approximate.  Neverdieless  the  method  gave  images  of  ~  1  cm  diameter  tumors 
in  roughly  60-70  %  of  the  cases  [Fantini  et  al.,  1996]  which  were  congruent  with 
X-ray  mammography  of  the  same  breast.  There  were  also  important  exceptions, 
but  the  results  are  the  first,  and  clearly  indicate  the  technique  is  potentially 
very  important  and  that  further  investigation  is  required. 

In  order  to  fully  realize  the  potentialities  of  the  optical  method  for  breast 
cancer  imaging  it  is  crucial  to  separate  spectroscopic  and  scattering  variations 
with  high  fidelity.  Imperfect  quantitation  will  severely  limit  both  the  sensitivity 
and  specificity  of  the  technique.  Complete  information  about  breast  tumor 
shapes  and  optical  properties  require  tomographic  approaches  to  the  data 
aquisition  and  inversion.  There  has  been  extensive  image  reconstruction  work 
with  diffusive  waves  along  these  lines  during  the  last  few  years  [see  for  example 
van  Houten  et  al.,  1996;  Ishii  et  al.,  1995;  Chang  et  al.,  1995;  Arridge  et  al.,  1991; 
Pogue  et  al.,  1995;  Oleary  et  al.,  1995].  In  essentially  all  of  these  approaches  a 
perturbative  model  (either  Bom  or  Rytov  approximation)  is  employed  to  define 
the  forward  problem.  Then  the  sample  volume  is  discretized,  voxel  weights  are 
assigned,  and  variations  in  absorption  and  scattering  factors  within  the  sample 
are  obtained  by  some  iteration  or  single- value  decomposition  procedure.  In  our 
lab,  for  example,  we  have  experimentally  demonstrated  the  tomographic 
reconstraction  of  absorptive  and  scattering  heterogenieties  in  tissue  phantoms 
[Oleary  et  al.,  1995;  O’Leary  SPIE],  we  have  shown  how  to  reconstract  lifetime 
and  concentration  of  fluorophores  in  tissues  [Oleary  et  al.,  1996a],  we  have 
reconstracted  the  dynamical  flow  properties  within  turbid  media  using  diffused 
temporal  light  field  correlation  functions  [Boas  et  al.,  1995],  and  we  have 
demonstrated  a  technique  whereby  stmctural  information  from  a  second 
imaging  modality  (such  as  X-Ray,  ultrasound,  or  MRI)  may  be  optimally 
combined  with  the  photon  migration  techniques  to  yield  more  accurate  optical 
properties  [Oleary  etal.,  1996b]. 

In  this  research  program,  we  proposed  to  develop  and  test  new  theoretical 
algorithms  for  breast  imaging  instruments  based  on  the  use  of  diffuse  photon 
density  waves.  One  scheme,  the  near-field  Fast  Fourier  Transform  (FFT) 
approach,  affords  a  rapid  method  for  imaging  of  tissue  interiors,  and  thus  will  be 
useful  for  on-line  analysis  of  the  breast  images.  The  second  scheme,  the  non- 
perturbative  approach,  is  a  novel  method  for  the  inclusion  of  boundary  effects, 
and  the  computation  of  absorption  and  scattering  factors  to  high  order; 
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importantly  this  scheme  does  not  involve  matrix  inversion  and,  as  a  result  of  its 
repeated  integrations,  is  expected  to  be  less  susceptible  to  experimental  noise.  If 
successful  these  powerful  theoretical  methods  are  expected  to  improve  the 
fidelity  of  the  breast  images,  and  it  is  anticipated  that  the  anatomical 
information  generated  by  the  device  will  be  useful  for  improving  our  ability  to 
specify  the  nature  of  breast  tumors. 


BODY 

We  have  implemented  several  aspects  of  the  proposed  work  during  the 
past  year.  Most  of  this  work  has  been  connected  with  the  near-field  FFT 
theoretical  approach,  although  all  of  the  instrument  development  will  be  useful 
to  both  theoretical  approaches.  We  briefly  review  these  developments  here. 

Experimental  Demonstration  of  Main  Theoretical  Ideas  Underlying  the  Near-field 
East  Fourier  Transform  Approach  rLi97.  Li981: 

In  order  to  demonstrate  the  feasibility  of  this  algorithm,  we  have 
performed  amplitde  and  phase  measurements  in  a  plane  parallel  geometry 
shown  schematically  in  Figure  1  below. 


Figure  1  : ;  Illustration  of  the  experimental  geometry:  the  source  is  fixed  at  the  origin 
and  the  detector  scans  in  a  planar  geometry  at  z=Zrf=5.0  cm  over  a  9.3  x  9.3  cm^ 
region  in  x-y  step  sizes  of  0.3  cm.  The  optical  heterogeneities  are  embedded  in  the 
highly  scattering  medium  between  the  source  and  the  detection  plane. 
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The  entire  experiment  was  embedded  in  an  effectively  infinite  sample  of 
Intralipid  so  that  infinite  media  boimdary  conditions  could  be  safely  employed. 

In  order  to  carry  out  these  experiments  we  developed  a  rapid  homodyne 
detection  systems  based  upon  In-phase/In-quadrature  demodulation 
techniques.  Much  of  the  experimental  and  theoretical  details  are  described  in 
our  publication  [Li97]  which  is  included  as  an  Appendix  to  this  report. 

The  demonstration  was  carried  out  using  ~3mW  of  light  at  786  nm.  Light 
from  the  diode  laser  was  amplitude  modulated  at  100  Mhz  and  fiber  coupled 
into  the  fishtank  containing  &e  Intralipid  (with  scattering  factor  of  about  7.3  cm’ 
*).  Scanning  of  the  detector  fiber  was  made  in  steps  of  about  0.3  cm.  The 
amplitude  and  phase  of  the  transmitted  diffuse  photon  density  waves  were  thus 
recorded  at  1024  points  on  the  detection  plane.  In  these  test  experiments  we 
first  measured  the  signals  with  a  uniform  background,  and  then  we  placed  two 
heterogeneities  (absorbing  or  scattering)  into  the  medium,  remade  the 
measurements,  and  extracted  the  projection  images  of  medium  optical  properties 
with  the  FFT  technique. 

Theoretically,  the  technique  consists  of  several  steps.  After  measuring  the 
“total”  diffuse  photon  density  waves  (due  to  the  heterogeneities  plus 
background),  we  subtracted  off  the  “background”  wave  to  obtain  the  “scattered 
wave”.  The  scattered  wave  is  then  spatially  fourier  transformed  in  the  plane  of 
detection.  The  FFT  method  formally  relates  the  fourier  transform  of  the 
scattered  wave  to  a  weighted  sum  of  the  fourier  transforms  of  the  so-called 
tumor  function;  the  tumor  function  essentially  contains  all  of  the  information 
about  the  sample  heterogeneities  such  as  their  position,  size  and  optical 
properties.  For  relatively  thin  heterogeneities  we  select  a  plane  (or  slice), 
parallel  to  the  detection  plane  but  within  the  sample  volume,  that  best 
approximates  the  position  of  the  heterogeneities.  Then  a  simple  algebraic 
operation  and  an  inverse-FFT  generates  an  two-dimensional  projection  image. 

We  demonstrated  this  method  in  [Li97],  and  the  results  for  two  absorbing 
heterogeneities  are  shown  in  Figure  2.  The  method  was  shown  to  work  for 
absorbing  and  scattering  objects  [Li97].  However  for  extended  objects  (i.e.  with 
diameters  larger  than  a  few  random  walk  steps)  we  found  that  the  absolute 
optical  properties  were  incorrect,  though  the  relative  optical  properties  were  still 
quite  accurate.  This  is  due  to  the  breakdown  of  the  thin  slice  approximation,  and 
requires  three-dimensional  reconstruction  techniques. 
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(C) 


Figure  2  I  (a)  shown  the  exact  x-y  positions  of  the  two  thin  absorbing  slices,  (b) 
shows  the  reconstructed  2-D  photographic  image  of  these  two  slices  using  the  K-space 
spectrum  analysis  algorithm,  (c)  shows  the  surface  plots  of  the  reconstructed  images 
where  the  reconstructed  absorption  variations  can  be  found. 
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New  Instrumentation  for  More  Realistic  Compressed  Breast  Geometries: 

There  are  many  improvements  that  must  be  made  upon  the  work 
described  above  [Li97]  before  the  FFT  methods  can  be  moved  to  the  clinic.  First 
and  foremost,  we  need  an  instrument  designed  to  work  in  a  realistic  geometry. 
Secondly  we  need  more  rapid  and  flexible  scanning  at  multiple  optical 
wavelengths. 

Over  the  past  six  months  we  have  developed  such  a  device.  We  designed 
the  instrument  onto  a  bed-like  frame  built  of  light-weight  aluminum  strutting. 
The  table  is  lockable  and  also  has  removable  wheels.  The  rigid  frame  is  critical 
to  eliminate  positioning  errors.  We  have  built  into  the  frame  a  two-plate,  soft- 
compression  optical  mammogram  unit.  One  of  the  plates  has  a  source  grid  of  12 
fibers  embedded  into  the  compression  plate.  By  comparison  to  our  single  source 
systems,  we  have  more  flexibility  to  define  our  scan  region,  and  we  anticipate 
that  we  will  need  multiple  source  positions  for  the  true  three-dimensional 
reconstructions.  The  other  critical  problem  we  have  addressed  in  the  current 
instrument  is  offset  drift.  In  particular  for  every  measurement  we  optically 
multiplex  between  a  750  nm,  an  786  nm,  and  a  “dark”  measurement.  This 
allows  us  to  continually  subtract  off  electronic  DC  offsets. 

On  the  detection  side  we  still  use  stepper  motors  to  scan  a  9cm  x  9cm 
detection  plate  in  0.6  cm  increments  (flexible).  The  previous  scanning  system 
was  fragile  at  best.  The  current  one  is  much  more  robust  In  addition  we  have 
added  in  an  electronically  controllable  attenuator  into  the  signal  path  before  the 
demodulator  in  order  to  increase  the  linear  dynamic  range  of  the  device. 

In  addition  to  building  the  new  instrument  we  have  carried  out  a  series  of 
meticulous  tests  in  order  to  correct  for  systematic  amplitude  and  phase  shifts  in 
the  electronics  and  to  discover  our  limitations  due  to  positional  noise.  These 
phantom  tests  are  underway  in  parallel  with  theoretical  improvements. 

Theoretical  Studies: 

We  have  implemented  several  improvements  in  our  theoretical  approach. 
The  first  problem  we  have  addressed  concerns  the  true  slab  geometry.  In 
particular  we  have  determined  and  incorporated  the  slab  geometry  Green’s 
Functions  into  the  original  FFT  theory.  This  has  lead,  for  example,  to 
considerably  better  agreement  with  analytical  theory.  The  detailed  theoretical 
formulas  are  extensive,  but  may  be  found  in  Xingdi  Li’s  thesis  [Li98]. 

Another  problem  that  is  often  glossed  over  in  the  literature,  but  that  we 
have  found  to  be  very  important  is  the  use  of  processing  “filters.”  We  have 
initiated  a  systematic  investigation  of  the  effects  of  various  filters.  Generally  we 
find  that  the  use  of  some  sort  of  filtering  is  crucial  to  reduce  reconstruction 
artifacts.  Examples  of  filters  used  are:  low  pass  filters  on  the  detection  (k- 
space)  grid  data,  Hamming  filters  on  the  tumor  function  (perturbation),  and 
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Hamming  filters  on  the  reconstructed  images  in  real  space.  Unfortunately  each 
of  these  filters  also  introduces  artifacts.  For  example,  the  reconstructed  optical 
properties  are  generally  reduced  by  50%.  Our  hope  is  to  understand  these 
artifacts,  correct  for  them,  and  thus  optimally  use  the  processing  filters. 

The  last  critical  theoretical  issue  concerns  full  three-dimensional 
reconstruction  with  the  methods.  This  requires  the  greatest  amount  of  new 
insight.  We  have  formally  developed  (unpublished)  an  approach  which  uses 
multiple  input  plane  waves  to  generate  3-D  images.  Of  course  plane  wave 
sources  are  not  very  realistic.  On  the  other  hand  we  believe  the  our  new  12 
source  compression  plate  offers  the  possibility  for  a  related  3-D  approach, 
because  each  point  source  can  be  considered  to  be  built  from  many  plane  waves. 
Much  theoretical  work  remains  on  this  for  the  next  year. 

CONCLUSIONS 

We  have  accomplished  a  lot  since  submission  of  this  proposal.  First  we 
demonstrated  by  experiment  [Li97]  the  validity  of  the  FFT  image  reconstruction 
approach.  These  experiments  were  carried  out  under  ideal  (i.e.  infinite  media, 
etc.)  conditions  however.  To  this  end  we  have  put  together  that  works  in  more 
realistic  scenerios  such  as  the  compressed  breast  geometry.  This  apparatus  will 
be  useful  for  realistic  phantom  studies,  and  eventually  (not  in  this  granting 
period)  for  optical  mammography  in  the  clinic.  The  new  apparatus  has  solved 
many  of  the  early  technical  problems  we  had  with  the  old  instrument  including 
the  problem  of  IQ-chip  offset  drift,  multiple  wavelengths  for  spectroscopy  (we 
now  have  multiple  wavelengths),  and  many  mechanical  instabilities  are  now 
corrected.  On  the  theoretical  side  we  have  implemented  slab  Green’s  Functions 
into  the  FFT  treatment,  and  are  systematically  exploring  the  use  of  processing 
filters.  We  are  just  beginning  to  tackle  the  full  3-D  FFT  theory;  this  will  be  our 
most  difficult  task  for  the  next  funding  period.  No  work  was  done  on  the 
Auxilliary  function  approach,  becuase  the  other  project  was  going  so  well.  We 
hope  to  consider  this  theory  also  in  the  next  year. 
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Recently  the  application  of  near-IR  diffusing  light  for 
biomedical  diagnosis  and  imaging  has  gained  favor  be¬ 
cause  of  its  noninvasive  nature,  economy,  and  novel 
contrast  relative  to  other  diagnostics.^"^  To  this  end 
a  variety  of  techniques  for  imaging  with  diffuse  light 
have  been  explored.^’^  Most  methods  use  direct  ma¬ 
trix  inversion  (e.g.,  singular-value  decomposition)  or 
iterative  techniques  (e.g.,  algebraic  reconstruction- 
simultaneous  iterative  reconstruction  technique,®  con¬ 
jugate  gradient  descent)  for  image  reconstruction. 

Here  we  introduce  a  new  near-field,  diffusive-wave¬ 
imaging  methodology,  using  techniques  similar  to 
those  of  conventional  diffraction  tomography.^  The 
near-field  diffraction  tomography  method  has  at¬ 
tracted  the  attention  of  a  few  researchers  in  the 
photon-migration  field.®  In  this  Letter  present  a 
rigorous  account  of  the  theory  and  provide  the  first 
experimental  images  of  absorbing  and  scattering  ob¬ 
jects  in  turbid  media  obtained  by  this  approach.  The 
method  differs  from  least-squares  techniques®  in  that 
it  is  fast  and  noniterative.  In  addition  to  providing 
information  about  the  position  and  shape  of  a  hidden 
object  or  objects,  pro  jection  images  can  be  used  to  de¬ 
duce  the  optical  properties  of  heterogeneities  without 
the  need  for  complex  reconstruction  procedures  such 
as  matrix  inversion  when  the  heterogeneities  are  thin 
and  information  on  their  depth  is  available. 

We  adopt  the  frequency-domain  picture  for  our  dis¬ 
cussion.  An  intensity  sinusoidally  modulated  light 
source  coupled  into  a  highly  scattering  medium  such  as 
tissue  produces  a  diffiise-photon  density  wave  (DDPW) 
that  propagates  outward  from  the  source.®  The  am¬ 
plitude  and  phase  of  this  DPDW  depend  on  the  absorp¬ 
tion  and  the  scattering  coefficients  within  the  turbid 
medium. 

^  In  uniform  turbid  media  the  DPDW  from  a  point 

source  at  Vs  detected  at  position  ra  has  the  form 
?^o(i'd,rs,6^)  =  Moexpiiko\r  -  rsl)/(47rDo|r  -  rgi). 
Here  co  is  the  angular  source-modulation  frequency. 
Mo  is  the  ac  amplitude  of  the  source.  Do  =  v/iSpso') 
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is  the  diffusion  coefficient  (t;  is  the  speed  of  light  in 
the  medium),  ko  =  [(— is  the  DPDW 
wave  number,  and  p^o  and  pso'  are,  respectively,  the 
homogeneous  absorption  and  the  reduced  scattering 
coefficients  of  the  medium.  In  heterogeneous  media, 
the  total  DPDW,  C7^(rd,  rs,  w),  is  a  superposition  of 
incident  [[/oC^d,  and  scattered  [[/i(rd,rs,  cy)] 

DPDW’s.  To  the  first  order  in  the  variation  of  optical 
absorption  and  reduced  scattering  coefficients,  the 
scattered  wave  is 

Ui(ra,rs,co)=  f  T[Uo(r,r,,co)]G(lra  -  vl)d^r .  (1) 

Jv 

Here  G(|rd  -  r|)  is  the  Green's  function  for  the  DPDW 
in  the  homogeneous  medium.  We  call  T[r,  C/q 
(r,rs,<y)]  the  inhomogeneity  function.  Tabs[r,  t/o 
(r,rs,<y)]  =  (Spa(r)v/Do)Uo(r,rs,co)  for  absorbing 
objects,  and  Tscatt[r,  G'o(r,rs,  &))]  ^  [Sps'(r)3Doko^/ 
v)  -  V  IniSps'  +  MsoO]  •  VJ7o(r,rs,<y)  for  scattering 
objects.  The  integral  is  over  the  sample  volume  V, 

We  consider  a  parallel-plane  geometry  [Fig.  1(a)] 
that  is  potentially  applicable  to  the  compressed-breast 
configuration.  For  this  case,  a  natural  basis  set  for  the 
Green’s  function  in  Eq.  (1)  is  the  simple  Weyl  expansion 
form  in  terms  of  spatial  frequencies  p,  q,  m}^: 

G(|r.  -  r|)  -  ^  ^  .xptipta  -  «) 

+  iqiyd  -  y)  +  im{zd  -  zf] ,  (2) 

where  we  assume  that  Zd>  z  without  losing  generality; 
m  =  {ko^  -  p2  _  ^2)1/2  >  0.  Using  Eq.  (2) 

and  taking  the  two-dimensional  (2D)  spatial  Fourier 
transform  of  both  sides  of  Eq.  (1)  [with  respect  to 
treinsverse  (x,y)  coordinates],  we  obtain 

Ui{p,q,Zd,rs,o))  =  —  J  T(p,q,z,rs,o)) 

X  exp\im{zd  -  zfldz .  (3) 

The  left  side  of  Eq.  (3)  is  the  2D  Fourier  transform  of 
©  1997  Optical  Society  of  America 
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(b) 

Fig,  1.  (a)  Experimental  geometry:  The  source  is  fixed 
at  the  origin,  and  the  detector  scans  in  a  planar  geometry 
at  z  =  Zd  ^  5.0  cm  over  a  9.3  cm  X  9.3  cm  region,  (b) 
Setup:  APD,  avalanche  photodiode;  A,  amplifier;  SSB I/Q, 
single-sideband  in-phase/quadrature-phase  demodulator; 
LO,  local  port;  RF,  radiofrequency  port;  ADC,  analog- 
digital  converter;  GPIB,  general  purpose  interface  bus. 

the  scattered  DPDW  measured  on  the  detection  plane 
z  =  Zd-  This  integral  equation  is  approximated  by 
N  ^  ' 

X  T(p,q,zj,rs,0))exp{-mzj) 

■'■=1  2m 

=  e^p(-imzd)Ui{p,q,Zd,rs,o)),  (4) 

where  Sz  is  the  discretized  step  size  and  N  is  the  total 
number  of  slices  in  the  z  direction. 

For  the  projection  image,  we  replace  zj  on  the  left 
side  of  Eq.  (4)  with  the  estimated  slice  position  of  the 
object.  We  drop  the  sum  over  all  other  zfs  and  then 
perform  a  2D  inverse  Fourier  transform  of  f  to  ob¬ 
tain  the  projection  image.  When  the  object  thickness 
is  of  the  order  of  several  transport  mean  free  paths 
[l/(/>ta  +  we  can  deduce  accurately  the  optical 

properties  of  the  object  or  objects.  For  thicker  objects 
(i.e.,  >5  mm),  the  average  over  the  size  of  the  object 
weighted  by  the  sum  of  exponential  phase  factors  re¬ 
duces  the  accuracy  of  the  optical  properties.  However, 
position  information  is  accurate,  and  the  relative  opti¬ 
cal  properties  of  multiple  objects  are  also  very  accurate. 

To  demonstrate  the  feasibility  of  this  algorithm,  we 
performed  amplitude  and  phase  measurements  in  a 
parallel-plane  geometry  [Fig.  1(a)].  The  experimental 
setup  is  shown  in  Fig.  1(b).  The  system  consists  of 
a  rf-modulated  (100-MHz),  low-power  ('-'3-mW)  diode 
laser  operating  at  786  mri.  The  source  light  is  fiber 
guided  into  a  large  fish  tank  of  50L  0.75%  Intralipid 
{pa  ^  0.02  cm“^,  =  8.0  cm“^),  enabling  us  to 

use  infinite  medium  boundary  conditions.  A  detection 
fiber  couples  the  detected  diffusive  wave  to  a  fast 


avalanche  photodiode.  A  single-sideband  demodulator 
is  used  to  find  the  homod3me  of  the  signal  and  the 
reference  wave  at  100  MHz.  The  dynamic  range  of 
our  current  setup  is  -2500.  The  source  and  detection 
fiber  optics  are  moved  by  automated  stepper  motors. 
The  system  is  very  stable  for  —3  h. 

The  experimental  geometry  is  shown  in  Fig.  1(a). 
The  source  position  was  fixed  and  taken  to  be  the 
origin  of  our  coordinate  system.  As  shown  in  Fig.  1(a), 
we  "‘made”  the  detection  plane  by  scanning  a  single 
detection  fiber  over  a  square  region  from  (-4.65, 
-4.65,  5.0)  cm  to  (4.65,  4.65,  5.0)  cm  in  a  plane  at 
Zd  =  5.0  cm  in  steps  of  size  Ax  =  Ay  =  0.3  cm.  The 
amplitude  and  phase  of  the  DPDW  were  recorded  at 
each  position  for  a  total  of  1024  points.  We  directly 
measured  the  amplitude  and  phase  in  the  homogeneous 
medium  to  obtain  C/o(*’dj  **8,  (o). 

Two  absorbing  slices,  each  1.5  cm  X  1.5  cm  X  0.4  cm 
were  then  submerged  in  the  turbid  medium  (0.75% 
Intralipid).  The  slices  were  made  of  resin  plus  Ti02 
and  absorbing  dye.  Slice  1  with  pai  =  0.20  cm”^  was 
placed  at  position  (’’1.6,  -0.3, 3.0)  cm,  and  slice  2  with 
Pa2  =  0.10  cm~^  was  placed  at  (1.6,  0.3,  3.0)  cm.  The 
scattering  coefficients  of  these  two  slices  are  the  same 
as  that  of  the  background,  i.e.,  8.0  cm“^.  We  obtained 
the  scattered  wave  co)  by  subtracting  the  ho¬ 

mogeneous  background  DPDW  i7o(rd,rs,  w)  from  the 
measured  signal  [//(rd,  Vs,  (o).  The  2D  Fourier  trans¬ 
form  t7i(rd,rs,(y)  leads  to  the  inhomogeneity  function 
T{p,q,Zj,rsi<o)  in  Eq.  (4).  We  use  a  priori  informa¬ 
tion  about  the  object  position(s)  in  the  z  direction  to 
select  a  single  image  slice,  e.g.,  a  slice  at  ^  =  2:0b j, 
where  the  inhomogeneity  function  is  T(  p,  q,  2obj,  rg,  d>). 
The  2D  inverse  Fourier  transform  of  T(  p,  q,  Zobj,  rg,  (o) 
gives  an  accurate  spatial  map  of  the  absorption  varia¬ 
tions  5p'a(jC,y,2obj)- 

The  reconstructed  images  from  the  experimental 
data  are  shown  in  Fig.  2.  The  two  objects  are  well  re¬ 
solved  (with  a  peak-trough  of  >2).  The  reconstructed 
x-y  positions  of  these  two  slices  are  approximately 
(—1.80, —0.25)  cm  and  (1.85, 0.25)  cm,  close  to  their 
true  x-y  positions  of  (-1.6,  -0.3)  cm  and  (1.6, 0,3)  cm. 
Inaccuracies  in  the  position  measurements  might  ac¬ 
count  for  the  discrepancy.  Images  shown  in  Figs.  2(b) 
and  2(d)  are  unprocessed.  The  reconstructed  absorp¬ 
tion  coefficients  are  well  above  the  background-noise 
level  and  close  up  to  the  true  values  [e.g.,  M^!^-(0.22± 
0.03)  cm-i  and  «  (0.13  ±  0.03)  cm“M.  The  noise 
in  the  reconstruction  is  mainly  from  the  finite  step  size 
and  scan  region,  the  positional  error,  and  the  electron¬ 
ics.  The  refractive-index  mismatch  between  the  ob¬ 
ject  and  the  background  medium  also  contributes  to 
the  inaccuracy  in  the  reconstructed  optical  properties. 
The  fast  Fourier  transform  (FFT)  calculation  takes  less 
than  200  ms  of  CPU  time  on  a  SunSparclO  worksta¬ 
tion,  and  the  presence  of  multiple  objects  does  not  in¬ 
crease  the  computation  complexity. 

The  feasibility  of  the  FFT  algorithm  for  imaging 
a  scattering  object  was  also  demonstrated.  The  ge¬ 
ometry  is  the  same  as  in  Fig.  1(a).  In  this  case 
we  used  spherical  objects  as  heterogeneities  instead 
of  thin  shces  to  test  the  accuracy  of  our  algorithm 
for  imaging  extended  objects.  Sphere  1,  of  radius 
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Fig.  2.  (a)  Slice  showing  the  exact  posi¬ 
tions  of  the  two  slices;  for  experiments 
with  scattering  objects,  two  spheres  are 
placed  so  that  their  centers  coincide  with 
the  slice  centers,  (b)  Projection  image  of 
slices,  (c)  Projection  image  of  spheres. 
The  recovered  absorption  (d)  and  scatter¬ 
ing  (e)  properties  of  the  slices  and  spheres, 
respectively,  are  also  shown. 


0.75  cm  with  «  16  was  placed  at  (-1.6, 

*“0.3, 3.0)  cm,  sphere  2,  of  the  same  radius  with 
^  26  cm"\  was  placed  at  (1.6,  0.3,  3.0)  cm. 
Both  spheres  have  the  same  absorption  coefficient 
as  the  background,  i.e.,  0.02  cm“^  We  see  from 
Fig.  2(c)  and  2(e)  that  positional  information  about 
these  two  scattering  objects  is  recovered.  In  this 
case  we  use  a  slice  through  the  sphere  center  at 
z  =  Zohj  and  obtain  a  2D  scattering  contrast  im- 
age.  We  do  not  expect  to  reconstruct  the  scattering 
wv  coefficients  accurately  since  the  objects  are  extended. 
However,  we  were  still  able  to  obtain  the  correct  posi¬ 
tions,  and  the  two  objects  are  well  resolved  [Figs.  2(c) 
and  2(e)] .  Furthermore,  the  ratio  of  the  reconstructed 
scattering  coefficients  is  close  to  the  true  ratio,  i.e., 
f^s2  / f^si  ^  1.35,  while  the  true  ratio  is  —1.62. 

Depth  information  is  required  for  full  3D  images 
with  this  diffraction  tomography  technique.  One 
simple  method  is  to  use  a  secondary  localization 
scheme  to  deduce  object  depth.  Alternatively,  two 
projection  images  of  the  sample  along  orthogonal 
directions  provide  sufficient  information  for  3D 
reconstruction. 


a 


We  have  successfully  applied  near-field  diffraction 
tomography,  diffuse-photon  density  waves,  and  FFT’s 
to  obtain  projection  images  of  hidden  objects  in  highly 
scattering  tissue  phantoms.  It  may  be  possible  to  ob¬ 
tain  clinical  projection  images  in  real  time  with  this 
^FT  approach.  The  geometry  used  so  far  has  been 
infinite.  In  practice,  bovmdary  effects  present  impor¬ 
tant  problems.  On  the  one  hand,  matching  materials 
might  be  used  to  reduce  the  boundary  effects.  On  the 
other  hand,  the  introduction  of  a  surface  integral  term 
in  Eq.  (1)  or  better  Green’s  functions  (which  vanish  on 
the  extrapolated  boundary)  may  be  used  to  incorporate 
boundary  effects.  The  technique  presented  here  pro¬ 
vides  a  basis  for  more  complicated  and  realistic  recon¬ 
struction  methods  to  address  these  issues. 
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